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We report on progress in our work on obtaining high precision standard model predictions for the Higgsstrahlung 
reaction at a linear collider. 



1. MOTIVATION 

Main production mechanisms of the standard 
model (SM) Higgs boson in e'^e~ collisions at a 
linear collider [1] are the Higgsstrahlung reaction 

e+e" ^ ZH (1) 

and the WW fusion process 

e+e- W*W* Ve>^eH. (2) 

The cross section of (1) decreases as 1/s while 
that of (2) grows as ln(s/m|^). Hence, while the 
Higgsstrahlung dominates the Higgs boson pro- 
duction at low energies, the WW fusion process 
overtakes it at higher energies. For small val- 
ues of the Higgs boson mass, vhh < 140 GeV, 
as is the case considered in this talk, the Higgs 
boson dominantly decays into a bb quark pair. 
As the Z boson of reaction (1) decays into a 
fermion-antifermion pair too, one actually ob- 
serves the Higgsstrahlung through reactions with 
four fermions in the final state. 

In the following we will concentrate on one spe- 
cific four fermion channel relevant for detection of 

(1) 

e+(pi) + e^{p2) M^ba) + M"(P4) 

+ Hp5) + bip6), (3) 
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where the particle four momenta have been in- 
dicated in parentheses, and the corresponding 
bremsstrahlung reaction 

e^e^ fi'^fi^bbj. (4) 

Reactions (3) and (4) should leave a particularly 
clear signature in a detector. To the lowest order 
of the SM, in the unitary gauge and neglecting the 
Higgs boson coupling to the electron, reactions 
(3) and (4) receive contributions from 34 and 236 
Feynman diagrams, respectively. Typical exam- 
ples of the Feynman diagrams of reaction (3) are 
depicted in Fig. 1. The Higgsstrahlung 'signal' 
diagram is shown in Fig. la, while the diagrams 
in Figs, lb and Ic represent typical 'background' 
diagrams. The Feynman diagrams of reaction (4) 
are obtained from those of reaction (3) by attach- 
ing an external photon line to each electrically 
charged particle line. 




Figure 1. Examples of the Feynman diagrams of 
reaction (3): (a) the double resonance 'signal', (b) 
and (c) the 'background' diagrams. 



In order to match the precision of data from 
a high luminosity linear collider, expected to be 
better than 1%, it is necessary to include radia- 
tive corrections in the SM predictions. At the 
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moment, however, it is not feasible to include the 
complete 0{a) corrections to a four-fermion re- 
action like (3) into a Monte Carlo (MC) gener- 
ator. Therefore it seems to be natural to apply 
the double pole approximation (DPA), as it had 
been done before in literature in case of the W- 
pair production [2] at LEP2. 

To the lowest order of the SM, the cross sec- 
tions of reactions (3) and (4) can be computed 
with a program ee4fg [3]. On the basis of ee4fg, 
we have written a dedicated program eezh4f that 
eventually should include all the leading radia- 
tive corrections to the four-fermion reactions rel- 
evant for the Higgs boson production and decay 
through the Higgsstrahlung mechanism. How- 
ever, up to now it includes the electroweak 0{a) 
radiative corrections to the on-shell reaction (1) 
of [4] and the initial state radiation (ISR) QED 
corrections. No corrections to Z and Higgs decay 
widths, nor higher order ISR effects, nor any non- 
factorizable corrections are included in eezh4f 
yet. 

2. OUTLINE OF THE CALCULATION 

The amplitude of (1) is represented in the fol- 
lowing way 

MzH {pi,o-i,p2,a2,kz,Xz,kH) = w (pi,(Ti) 
X {pi,P2, kz, kn) e* ikz,Xz) u (p2, cr2) , (5) 

where pi,P2,kz,kH and ai,a2,Xz are the on- 
shell four momenta and helicities of the external 
particles. Neglecting the electron mass, the am- 
plitude {pi,p2,kz,kH) is decomposed in the 
following Lorentz covariant way 

F^' ipi,p2, kz,kH) = AiY' + ^27^75 (6) 

+ (^3lPi + ^32P^) h + {AhpI + A^2p'i) 

where we have changed a little the original nota- 
tion with respect to that of Eq. (2.5) of [4]. The 
six invariant amplitudes that do not vanish in the 
limit me ^ are now denoted by Ai and Aji, 
i = 1.2, 7 = 3,4. In order to fix normalization 
we give here the two amplitudes, A^^\ i = 1,2, 
that are nonzero in the lowest order. They read 

(0) _ 9HZZ 



with the vector and axial-vector coupling of the 
Z to leptons 



ai = ew, 

'iSwCw 



a2 = 



ew 



4swCw 



(8) 



and the Z-Higgs boson coupling qhzz = 

ewmz /{swcw)] the electric charge ew = 
1/2 

(Anaw) and electroweak mixing angle 0w are 
fixed by 

aw = V2G^m^s'^/Tr, s^ = l-TO^/m|, (9) 

where mw and niz are the physical masses of 
the W and Z boson. This choice, which is ex- 
actly equivalent to the G^^-schemc of [4] , is in the 
program referred to to as the 'fixed width scheme ' 
(FWS). We have introduced a shorthand notation 
sw = sin 6w and Cw = cos 9w in Eqs. (8-9). 

Amplitudes Ai, i = 1,2, of Eq. (6) are split 
into the Born part of Eq. (7) and the one 
loop corrections in the following way 

Ai = ^ + MiQED + SAi, i = l,2, (10) 

where (J^iQED denotes the QED infra-red (IR) 
divergent and SAi the IR finite part of the one- 
loop electroweak corrections. The IR divergent 
part 6AiQED is factored out and combined with 
the soft bremsstrahlung correction to (1) result- 
ing in a QED initial state radiation correction 
factor Cq^y) [4] that can be split into a univer- 
sal cut-off dependent and non-universal constant 
part 

Cqed Ecut) = C'qed (^5 Ecut) + C'qed (11) 
with 



Cqed {s,Ecut) — 



In ^ - 1 In 



3, s 
4 mi 



(12) 



and 

^non— univ. ^ / _ i \ 

^QED -^ly ^i- 



(13) 



A) 



s — + imz^z 



(7) 



1 /2 

The electric charge e = (47rao) in Eqs. (12 13) 
is given in terms of ao in the Thomson limit. 
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The one-loop amplitudes 



SA^ = 6Ai{s,cos0), Aji = Aji{s, cos 9) (14) 

are complex functions of s = {pi +^2)^ and cos^, 
9 being the Higgs boson production angle with 
respect to the initial positron beam in the cen- 
tre of mass system (CMS). They are computed 
numerically with eezh4f making use of FF 2.0, 
a package to evaluate one- loop integrals by G. J. 
van Oldenborgh [5]. 

As the computation of the one-loop electroweak 
amplitudes of Eq. (14) slows down the MC inte- 
gration substantially, a simple interpolation rou- 
tine has been written that samples the amplitudes 
at a few hundred values of cos 9 and then the am- 
plitudes for all intermediate values of cos 9 are 
obtained by a linear interpolation. This gives 
a tremendous gain in the speed of computation, 
while there is practically no difference between 
the results obtained with the interpolation rou- 
tine and without it. 

The one-loop virtual and soft bremsstrahlung 
corrected matrix element squared of (3) in the 
double pole approximation (DPA) reads 



ISR 



virt-t-soft 



(15) 



M 



(0) 

DPA 



r" 



non— univ. 



QED 



2Re (^M^°'>p^6Mdpa 



where M^^^ is the lowest order matrix element of 
(3) and the QED correction factors are given by 
Eqs. (12) and (13); the lowest order matrix ele- 
ment M^jpj^ and the one- loop correction 5Mupa 
in the double pole approximation for the Z and 
Higgs boson are given by 



M 



(0) 



DPA - vi'Pi) {^l^^f^ + ^2°^7/x75) u (pa) 

X P^Ph (16) 

SMdpA = V (pi) [5^17;:, + ^27^75 
+ (-431^1;:^ + ^32P2/x) h 
+ {AiiPif, + A42P2n) his] U (P2) 

X P| Ph. (17) 

The pole factors P| and Ph in Eqs. (16) and (17) 
are given by 



P^ 



u{ki) jari' +(12-''-:,) i:{kz) 



Ph 



9Hbbu{h)v{k5) 



{P5 + Pe) - rii\ + imnVH ' 



(18) 



with the Z-fi couplings ai , 02 given by Eq. (8) 
and the Higgs-bottom-quark coupling gnhb = 
-ew/{'^sw){rnb/mw)- 

The projected momenta kz and ki, of Eqs. (17) 
and (18) are obtained from the four momenta Pi, 
i = 3, 6, of the final state fermions of reaction 
(3) with the following projection procedure. 

First the on-shell momentum and energy of the 
Higgs and Z boson in the CMS are found 



Ai/2 



2si/2 



kn = 



kn 



P5+P6 

|P5 +P6| 



E 



H 



kz 



-kH, Ez = \/s-E, 



H- 



Denote four momenta P5 and ps of reaction (3) 
that are boosted to the rest frame of the Higgs 
and Z, respectively, by pg and P3. The projected 
four momenta fc^, i = 5, 6, in the rest frame of the 
Higgs are obtained with 



AV2 [r 



''Hi "'55 



5, 



2mH 



E'i 



kU = 



P'5 



P'5 



Similarly one obtains k'^ and k'^ in the rest frame 

of the Z. 

Four momenta k\, i = 3, 6, are then boosted 
back to the CMS giving the projected four mo- 
menta ki, i = 3, 6 of ^'^ ,11^^ ,h and b that sat- 
isfy the necessary on-shell relations 



fcg 



— ^4 



kl 



(fcs + fcg) = m 



H- 



The actual value of cos 9 in Eq. (14) is given by 
cos^ = k%/\kH\- The described projection pro- 
cedure is not unique. As the Higgs boson width 
is small, the ambiguity between different possible 
projections is mainly related to the off-shellness 
of the Z boson and is of the order of aTz/i'Kmz)- 

3. NUMERICAL RESULTS 

In this section, we present some numerical re- 
sults for (3) which will be one of the best detection 
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channels of the Higgsstrahlung reaction (1) at the 
Hnear coUider. 

The computation has been performed with a 
program eezh4f in the FWS with the following 
set of initial SM electroweak physical parameters 
[7]: 

mw = 80.423 GeV, Tw = 2.118 GeV, (19) 
mz = 91.1876 GeV, Tz = 2.4952 GeV, (20) 
= 1.16639 X 10"^ GeV"^ (21) 
ao = 1/137.03599976. (22) 

The charged lepton masses in MeV are 

me = 0.510998902, = 105.658357, 

rur = 1776.99 (23) 

and the heavy quark masses in GeV are 

rUc = 1.5, mt = 174.3, rub = 4.4. (24) 

The light quark masses in MeV are assumed at 

m„ = 62, rud = 83, = 215, (25) 

which together with as = 0.133 reproduces the 
hadronic contribution to the running of the fine 
structure constant. The Higgs boson mass is as- 
sumed to be niH = 115 GeV and its width is 
calculated to the lowest order of the SM. 
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Figure 2. The 'signal' total cross section of reac- 
tion (3) in the narrow width approximation as a 
function of CMS energy. 



In Fig. 2, we plot the total 'signal' cross sec- 
tion of (3) in the narrow width approximation, 
i.e. the cross section of (1) multiplied by the cor- 
responding branching ratios, for the Higgs and 
Z boson, as a function of the CMS energy. The 
dashed curve shows the Born cross section while 
the solid curve shows the cross section includ- 
ing the one-loop electroweak corrections without 
the QED initial state radiation correction factor 
^QED of Eq. (11). The corresponding relative cor- 
rection is plotted in Fig. 3. The correction is of 
the order of a few per cent for ^/s < 500 GeV and 
it grows almost linearly in the absolute value for 
higher energies, most probably due to large Su- 
dakov logs that enter through the wave function 
normalization factors. Once the corrections to 
the Z and Higgs boson decay are included, these 
largo corrections combined with the wave func- 
tion renormalization of the Z and Higgs decay 
width should cancel against the renormalization 
constant of the Z and Higgs propagator. 
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Figure 3. The relative corrections corresponding 
to that of Fig. 2 



In Fig. 4. we compare results for the total cross 
sections of (3) as a function of the CMS energy, 
obtained with a complete set of the Feynman di- 
agrams in the lowest order of the SM, including 
the Higgs boson (solid line) and without the Higgs 
boson exchange (dashed line). The Z-Higgs pro- 
duction signal is visible as the second bump of 
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the solid line, the first bump reflecting the dou- 
ble Z production resonance. The dotted line in 
Fig. 4 shows the 0(a) corrected cross section in- 
cluding the QED ISR correction (12) applied to 
the complete set of Feynman diagrams of (3), the 
initial state hard photon radiation integrated over 
the full photon phase space and the electroweak 
radiative corrections to (1) in DPA. Fig. 5 fo- 
cuses on the CMS energy region, where both res- 
onances are present. We see that the photon ra- 
diation smears the resonances, as expected. The 
somewhat artificial large logarithmic correction of 
Fig. 3 has practically no visible effect in Figs. 4 
and 5 as the 0{a) is dominated by the QED ISR 
correction. 
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Figure 4. Total cross section of (3) as a function 
of CMS energy. 



4. SUMMARY AND OUTLOOK 

The status of work on the high precision SM 
predictions for the Higgsstrahlung reaction at the 
linear collider has been reported. The QED ISR 
corrections to e+e~ ii^fi~bb have been in- 
cluded. The virtual one-loop weak corrections to 
(3) in the DPA have been included, too. Work 
on inclusion of the corrections to Z and Higgs 
decay widths and higher order corrections due to 
the QED ISR is in progress. The non-factorizable 
corrections should be included, too. 
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Figure 5. Total cross section of (3) in the vicinity 
of the ZZ and ZH resonances. 
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